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Abstract 

In this paper we analyse Abelian diagonal orthogonally transitive 
spacetimes with spacelike orbits for which the matter content is a 
stiff perfect fluid. The Einstein equations are cast in a suitable form 
for determining their geodesic completeness. A sufficient condition 
on the metric of these spacetimes is obtained, that is fairly easy to 
check and to implement in exact solutions. These results confirm that 
non-singular spacetimes are abundant among stiff fluid cosmologies. 

PACS 04.20.Dw, 04.20.Ex, 04.20.Jb 

I Introduction 

After the discovery of the first non-singular perfect fluid cosmological model 
by Senovilla ilj, the possibility of constructing regular cosmologies was re- 
newed. The interest for regular cosmologies had stifled for nearly thirty 
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years due to the powerful singularity theorems (cfr. for instance 0), 
which seemed to preclude such spacetimes under very general requirements, 
such as chronology protecting, energy and generic conditions. The open way 
to regular cosmologies was found in the violation of some technical premises 
of the theorems. For instance, in pi it was shown that the Senovilla space- 
time did not possess a compact achronal set without edge and could not have 
closed trapped surfaces. 

However, the first results were not encouraging. The extension of the 
Senovilla solution to a family of spacetimes left the set of regular models 
limited to a zero-measure subset surrounded by spacetimes with Ricci and 
Weyl curvature singularities 0. During the following decade only a few new 
non-singular cosmologies were added to the list 0. 

Another estrategy to approach singularities arose with the publication of 
regularity theorems [Zj, |H], [Oj- Whereas singularity theorems stated general 
sufficient conditions for the appearance of singularities, these theorems aimed 
the contrary, namely particular conditions to achieve regular spacetimes. 

The application of the conclusions of |H| to a restricted family of stiff 
fluids provided an unexpected result. The set of known non-singular perfect 
fluid cosmologies was enlarged with a huge family depending on two nearly 
arbitrary functions jlQj . 

The purpose of this paper is the extension of those results to determine 
which spacetimes among Abelian diagonal orthogonally transitive spacetimes 
with spacelike orbits and with a stiff fluid as matter content are non-singular. 
Instead of restricting to an integrable family of solutions of the Einstein equa- 
tions, we analyse the whole set of diagonal cylindrical stiff fluid spacetimes 
with a spacelike transitivity surface element. 

With this aim in mind we write in Section |n] the Einstein equations for 
such spacetimes and we cast them in a form suitable for the application of 
the theorems. The analysis of the restrictions imposed by regularity condi- 
tions is done in Section IIIIl Finally in Section IIVI we check the possibility 
of constructing regular spacetimes with non-vaninishing matter scalar space 
averages on Cauchy hypersurfaces in order to support the validity of a regu- 
larity conjecture by Senovilla [TTj . 

II Equations for G2 stiff fluid spacetimes 

As it has been stated in the introduction, we shall focus on spacetimes en- 
dowed with an Abelian orthogonally transitive group of isometries G2 acting 
on spacelike surfaces, since this is the framework where most non-singular 
spacetimes have been found so far. We further impose that generators for 
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the group can be found that are mutually orthogonal. We follow jHj for writ- 
ing the Einstein equations for such spacetimes using a formalism based on 
differential forms. 

If the generators of the isometry group are chosen to be {C,v}y we may 
write an orthonormal tetrad, {6^,6^,6^,6^}, where just 6"^ and 6^ lie in 
\m{^,ri}. We may impose that these 1-forms be Lie-invariant under the 
isometry group j^. The metric is written as 

ds'^ = -6^ (^e^ + 6^ ®6^ + 6"^ ®6'^ ^6^ ® 6^. (1) 

Making use of the spacelike congruence for 6"^ and its kinematical quan- 
tities, we may define the tetrad basis according to the vanishing torsion 
equations, 

c^^o = z/A^^ (2a) 

d6^ = uA6^ (2b) 

d6^ = a A 6'^ (2c) 

d6^ = {f3-a)A6\ (2d) 

where u is just a connection in the 6^ — 6^ subspace, a is an "acceleration" 
for 6*^ and (3 is related to the expansion of the surface element in the 6^ — 6^ 
subspace, since d{6'^ A6^) = (3 A6'^ A 6^. 

The integrability conditions for these equations are easily obtained by 
exterior differentiation of the system. 



dp = 0, (3a) 

da = 0. (3b) 

Finally, Einstein field equations are written in terms of these differential 
forms as an exterior system, 

d*a + pA*a=(^^T-T22^6^A6\ (4a) 

d*p + pA*P=iTu~Too)6^A6\ (4b) 

du + a A *a - f3 A *a = ^{Too - Tn + T22 + T33) 6^ A 6\ (4c) 

* /3 + /5 A */3 + 2(a - /3) A *a + 2z/ A /3 = (Too + Tn) 6^ A 6\ (4d) 

ci/3 + /3 A /3 + 2(a - /3) A a + 2i/ A */3 = 2Toi^° A 6^, (4e) 
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for a matter content defined by the energy-momentum tensor T = TabO"'®d^ . 

The tilde denotes a reflection in the 9^ — 9^ subspace, that is, if a = 
a9^ + 6^^, then a = a9^ — b9^. The * denotes the Hodge duality operator in 
the same subspace, *a = —a9^ — b9^. 

Integration of the first Bianchi equations (jHafHbj) . 



a = —dU, (5a) 

13 = , (5b) 

P 

allows integration of Cartan equations in terms of two functions, z, (p, 

that we take as coordinates in order to write the metric in a conventional 
form, 

g = e^^i-df + dr^) + e-^^'dz^ + p^e^^rf^l (7) 

The coordinates are adapted to the Killing fields, so that ^ = d^, r] = d^. 
The non-ignorable coordinates t, r are chosen so that the metric is isotropic 
in the subspace spanned by 6*° and 9^, 

00 = e^dt, 9^ = e^dr. (8) 

The range for these coordinates is the usual one, 

-cx) < t, 2; < cx), < r < cx), < < 27r, (9) 

if we require the spacetime to be cylindrically symmetric. The remaining 
metric functions, K, U and p, depend just on t and r. 
The connection in this case is u = *dK. 

This is the general framework for an orthogonally transitive diagonal 
spacetime with spacelike orbits. If the matter content is a perfect fluid with 
4-velocity m, pressure p and density /i, the Bianchi equations for such energy- 
momentum tensor, 

T = pu ® u + p {g + u ® u)^ (10) 
may be written in compact expressions involving the kynematical 1-forms, 

du-\ dpAu = 0, (11a) 

P + P 

d*u+ ( (3+ A*u = 0, (lib) 
V P + Pj 
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which state that the fluid is irrotationaL 

We might choose 9^ = u for writing the Einstein equations, as it was 
done in ^0], but since we aim full generality, we shall not follow that way 
and explore arbitrary possibilities of alignment for this 1-form. Preserving 
the unitarity of m, we may parametrize it in terms of a function ^, 

M = -0° cosh^-^^ sinh^, (12) 
so that the Einstein equations for a perfect fluid take the following form 

Utt - Urr + -{Utpt - UrPr) = ^^e^^, (13a) 

P 2 

Ptt - Prr = {p- P)pe'^^, (13b) 

KtPr + KrPt = Ptr + UtPr + UrPt + 2pUtUr 

+e2^p/i±^sinh2e, (13c) 



KtPt + KrPr = ^" ^ + UtPt + UrPr 

+p (^U^ + + e'^^ cosh2e) , (13d) 
- + ^^rP.^ + U^^ -U^ = ^ e-, (13e) 
and the energy-momentum conservation laws yield 

Pr cosh^ i+ {pt — Pt) sinh ^ cosh i — Pr sinh^ ^ 



P + P 



+ ^*^"^^^-^-^^"^^sinhe = 0, (14a) 

Pt cosh^ ^ + {pr — Pr) sinh ^ cosh ^ — pt sinh^ ^ 

J<t - H : 

P + P 

^p,coshe-p.sinhe^^^^ (14b) 

p 

The system of equations becomes much simpler if we restrict to stiff fluids. 
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p = p, 



Utt - Urr + -{UtPt - UrPr) = 0, (I5a) 
P 

Ptt - Prr =0, (15b) 
KtPr + Krpt ptr + UfPr + Uypt . ^jj jj . 2K • , oC 

= V2UtUr + e psinn24, (15c) 

P P 

KtPt + Krpr _ ptt + Prr UtPt + UrPr jj2 , jj2 

— + + + u 

p 2p p 

+e2^pcosh2^, (15d) 

I ^rPr ~ UtPt jj2 tt2 2K ^i c; N 

Krr - Ktt H VU,^ - U^ = pe , (15e) 

p 

J. c , P'- I Ptcosh^-p^sinh^ • c n ^^^(\ 

^r- - 6 + TT Sl^^ ^ = 0' ) 

2p p 

J. c , Pt , Pt cosh ^ - prSinh ^ c n \ 

Kt-^r + 7r + = 0. (15g) 

2p p 

The reason why the stiff ffuid equations are easy to integrate is that the 
metric functions U, p decouple from the pressure, which only appears in the 
equations for the conformal factor K. Therefore the stiff fluid case is fairly 
similar to vaccuum and can be generated from this one. 

A further simplification can be obtained if we take p as coordinate. This 
is fully compatible with an isotropic parametrization, since equation ()15b|l . 



^ ^ ^ , f *dp\ dpA*dp d * dp 
\ P J P P 

states that *dp is also an exact differential form. 

We take p = r as an spatial coordinate, since every known non-singular 
solution has a surface element with spacelike gradient. With this choice of 
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coordinates the differential system becomes even simpler, 



Utt- 


u 

Urr - — = 0, 

r 


(16a) 


K, — 


TL _)_ 2rUJJ 4- p^^rtr sinh 2f 




Kr = 


: Ur + r(f/2 + f/2) + e^^pr cosh 2^, 


(16c) 




r 


(16d) 


Kr - 


sinh^^ 
Zp r 


(16e) 


Kt- 


Pt sinh f cosh f 

+ ^ ^ ^ = 0. 

2p r 


(16f) 



The integrability condition, Krt = Ktr, for (jl6b[ I16cj) requires that a 
combination of functions be an exact differential form, 

dH = e^^rp{smh2^dt + cosh 2^ dr), (17) 
from which we can read ^ and the pressure, if K is known, 

TT p-2K I 

tanh2e = \p\ = ^Jh? - HI (18) 

The integrability of dH is also a consequence of the energy-momentum 
conservation equations ()16e[ I16fjl . 

For consistency these expressions imply that the gradient of H be space- 
like and that Hj. be positive in order to have positive pressure. 

The simple case, ^ = 0, for which u is parallel to the time direction 
corresponds to H = 7r^/2, where 7 is a positive constant. 

The remaining system of differential equations. 



Utt -Urr- — = 0, (19a) 
r 

r 

Kt = Ut + 2rUtUr + Ht, (19c) 

Kr = Ur + r{U^ + U^) + Hr, (19d) 

is formed by a reduced wave equation in polar coordinates for U on the plane 
and a non-linear wave equation for H. Once these equations are solved, we 
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are left with a quadrature for K. The integrabihty of this quadrature is 
guaranteed by the other equations. 

As it has aheady been stated, these equations are pretty similar to those 
of vaccuum. The only diference is the additional conformal factor defined by 
H. 

Regularity of the metric at the axis r = is already implicit in the 
equations, provided that metric functions are regular. Following [12], we 
have a regular axis whenever 

,.^^ (gradA gradA> ^^,(,_,,^^^^^^^ A = (9,. 9,) = r V. (20) 

r^O 4/\ 

But according to equations (fT71ll9cl|l9d|) . K and U are equal at the axis, 
except for a constant, since 

Kr{t, 0) = Ur{t, 0), Kt{t, 0) = Ut{t, 0), (21) 

if pressure and K are regular functions, so that Ht{t,0) = = iJ,.(t,0), and 
therefore condition (pUj) is fulfilled either by taking the constant equal to zero 
or conveniently rescaling the angular coordinate. 

The problem of obtaining solutions for H is solved by the Wainright- 
Ince-Marshman formalism ^3]. Solutions to ()19bj) may be generated from 
solutions of the reduced wave equation on the plane with timelike gradient, 

(Ttt -arr- — = 0, aj - al > 0, (22) 
r 

by a quadrature identical to the one which defines K — U in the vacuum case. 



Ht = 2rat(Tr, (23a) 
Hr = r{a^ + a^^). (23b) 

The functions H generated by this mechanism have trivially a spacelike 
gradient and positive radial derivative. The fluid properties may be read 
directly from the generating function, 

tanh 2e = P = e-'^'ia^ - a'r)- (24) 

The function that generates the ^ = case is a = ^/^t. 
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Using this formalism, the remaining system of equations is formed by a 
quadrature and two reduced wave equations, 



Uu - Urr - — = 0, 



(25a) 



r 



(^tt - o-fT = 0, 



(25b) 

(25c) 
(25d) 



r 



Kt = Ut + 2rUtUr + 2r(7t(Tr, 

Kr = Ur + r{U^ + U^) + r{al + al). 



Ill Non-singular models 

We have obtained a fairly simple system of equations ()25all25c|) that will be 
useful for analysing the regularity of the solutions. Following |i2j we take 
causal geodesic completeness as our definition for regularity. 

Even if we have regular metric components, geodesic completeness of 
the spacetime is not guaranteed and we have to check explicitly that every 
timelike and lightlike geodesic in the spacetime can be extended to all values 
of the affine parameter, that is, in the parametrization for which the geodesic 
equations take the form 



in terms of the Christoffel symbols. 

Fortunately, results concerning causal geodesic completeness of diagonal 
Abelian orthogonally transitive spacetimes have already been obtained in jH] . 
The conclusions of that paper may be condensed in two theorems. We follow 
the simplified version of |iin|| . 

Theorem 1: A diagonal Abelian orthogonally transitive spacetime with space- 
like orbits endowed with a metric in the form ((Zj) with metric 
functions K,U,p, where p has a spacelike gradient, is future causally 
geodesically complete provided that along causal geodesies: 

1. For large values of t and increasing r. 



(a) {K -U -\np)r + {K -U - In p)t > 0, and either {K - U - 
lnp)r > or \{K-U-\np)r\ < {K-U-ln p)r + iK-U-\n p)t. 

(b) Kr + Kt> 0, and either Kr > or \Kr\ < K.,- + Kf 

(c) {K + U)r + {K + U)t > 0, and either {K + U)r>Q oi\{K + 

U)r\<{K + U)r + {K + U)t. 



A^^ _|_ 

•Xj I J- j •Xj •Xj ^ 



(26) 
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2. For large values of t, a constant b exists such that 
K{t,r)-U{t,r) ] 

2K{t,r) \ > - ln\t\ + b. 

K {t, r) + U{t,r)+ In p{t,r) J 



Theorem 2: A diagonal Abelian orthogonally transitive spacetime with space- 
like orbits endowed with a metric in the form ((Tj) with metric func- 
tions K, U, p, where p has a spacelike gradient, is past causally geodesi- 
cally complete provided that along causal geodesies: 

1. For small values of t and increasing r, 



(a) {K-U - \np)r -{K-U- \np)t > 0, and either {K - U - 
Inp), > or \{K-U-\np)r\ < {K-U-\n p)r-iK-U-\n p)t. 

(b) Kr- Kt> 0, and either iT^ > or \Kr\ < Kr - Kf 

(c) {K + U)r -{K + U)t > 0, and either {K + f/), > or + 

U)r\<{K + U)r-{K + U)t. 

2. For small values of t, a constant b exists such that 



Therefore now we just have to verify under which conditions these the- 
orems can be applied to stiff fluid spacetimes. Since the theorems do not 
make use of Einstein equations, it is expected that when we take them into 
account the conditions will not be so restrictive as they seem. 

We begin with future-pointing geodesies. The first part of the theorem is 
a set of conditions on the derivatives of the metric functions: 

1. (a) From ()19cll9dj) we obtain 



The sum of the derivatives of H is always positive, since Hr > \Ht\ 
in order to have positive pressure. In fact, this is the r{at + (TrY 
term in the Wainright-Ince-Marshman formalism. 

This expression is positive if either \Ut + Ur \ or Ht + Hr {\crt + ar\ 
in the Wainright-Ince-Marshman formalism) does not decrease as 
1/r or faster for large values of t and r. That is, we need either U 
or H to overcome the negative term. 

Under such conditions, the second part of the premise. 



K{t,r) - U{t,r) 
2K{t,r) 
K{t,r) + U{t,r) + In p{t,r) 




{K-U - \np)t + {K-U - \np)r = r{Ut + Urf + Ht + Hr 



1 



r 
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{K-U - \np)r = r{U^ + f/^) + if, - - > 0, 

r 

is also satisfied. 

(b) Once ()la|l is fulfilled, this condition, 

Kt + Kr = Ut + Ur+ r{Ut + Urf + Ht + Hr > 0, 

is trivial, since the only possible negative contribution would be 
that of Ut + Ur and this is counteracted by Ht + Hr if it decreases 
as 1/r or faster, or by r{Ut + UrY if it does not. 

Following a similar line of thought we also conclude that 

is positive for large values of t and r. 

(c) The last set of conditions on the derivatives, 

{K + U)t + {K + U)r = 2{Ut + Ur) + r{Ut + Urf + Ht + Hr > 0, 

Kr + Ur = 2Ur + rU^ + Hr > 0, 

is also a consequence of (fTaj) . 
Therefore the first part of the theorem is satisfied if 

r^-'\Ur + Ut\ ^ 

or I 74 (27) 

r'-'{Hr + Ht) J 

for large values of t and r. 

The conclusion for past-pointing geodesies is quite similar. We just 
have to change the sign of the time derivatives. 

r'-'\Ur-Ut\ ] 

or I 74 (28) 

r^~%Hr-Ht) J 

for large values of r and small values of t. 

For instance, these restrictions are trivial for the = case, since 
H = 7r^/2 does not decrease. 
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2. The dependence on the matter content of the spacetime may be re- 
moved from these conditions, since we may write 

K{t,r) = U{t,r)+ [ dr' {r'U^{t,r') +r'Uf{t,r') + Hr{t,r')) 

Jo 

= f/(t,0)+ / dr'Krit,r'), (29) 
Jo 

and according to (jl9dj) or (j25d|l Kr is a positive term if the first part 
of the theorem is satisfied. 

(a) The first condition is tautological since 

K{t,r)-U{t,r) = / dr' {r'U^{t,r')+r'Uf{t,r') + Hrit,r')) >0. 
Jo 

(b) For geodesies along the axis, this condition requires for large values 
of the time coordinate that 

K{t,0) = U{t,0) > ~^\n\t\+b, (30) 

and for general geodesies the only difference is the positive term 
in ()29|) . Therefore ()30|) is the only restriction for all geodesies. 

(c) The same restriction is achieved likewise when applied to the ex- 
pression K + U -|- In p. 

Therefore we are left with just three regularity conditions on the metric of 
an Abelian diagonal orthogonally transitive spacetime with spacelike orbits 
and with a stiff perfect fluid as matter content. We may summarize these 
results in two theorems: 

Theorem 3: A cylindrical spacetime with a stiff perfect fluid as matter 
content, endowed with a metric in the form (|7j) with metric functions 
K, U, p is future geodesically complete if the gradient of the surface 
element is spacelike and 

1. For large values of t, U(t, 0) > — - In \t\ + b. 

2. Either r'^-^\Ur + Ut\ or ^(-f^r + Ht) does not tend to zero for 
large values of t and r. 
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Theorem 4: A cylindrical spacetime with a stiff perfect fluid as matter con- 
tent, endowed with a metric in the form ((7j) with metric functions 
K, U, p is past geodesically complete if the gradient of the surface ele- 
ment is spacelike and 

1. For small values of t, U (t, 0) > — - In \t\ + b. 

2. Either r^~^\Ur — Ut\ or r^~^{Hr — Ht) does not tend to zero for 
small values of t and large values of r. 

For vaccuum spacetimes both theorems hold just dropping the conditions 
on the derivatives of H. 

The restrictions imposed by both theorems in order to have a non-singular 
spacetime are rather simple to implement, since U is just a solution of the 
wave equation and H is related to another one. We may state that regularity 
conditions are quite weak for stiff fluids, since it is very easy to provide 
solutions that fulfill such requirements. For instance: 

Corolary: A metric with arbitrary H and a function U which grows for 
large \t\ and for large r makes the spacetime geodesically complete. 

It is not difficult to derive such functions. The solutions to the reduced 
wave equation in the plane can be written as solutions of the initial value 
problem, 

Utt — Urr — = 0, 

U{0,r) = f{r), Ut{0,r)=g{r). (31) 
The solution to this problem can be written in closed form |151| . 



-1 pZTT PL 

U{x,y,t) = — d<p I dRR 
2vr Jo Jo 



+ 2^oil 'n/'''' — — ' ^''^ 

for initial data U{x,y,0) = f{x,y), Ut{x,y,0) = g{x,y), taking into account 
that / and g are to have circular symmetry. 

If we split U in Uf and Ug, the terms depending respectively on the 
initial data for U and its derivative, we notice that Uf is even in the time 
coordinate whereas Ug is odd. This implies that if Uf satisfies the first 
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condition in theorem 3 for future-pointing geodesies, it will fulfill it for past- 
pointing geodesies too. On the contrary, if Ug satisfies it for future-pointing 
geodesies, it may not fulfill it for past-pointing ones, unless it behaves for 
large values of |t| slower than a logarithm. Therefore one is to require either 
that Uf dominates over Ug for large |t| or that both terms behave slower than 
a logarithm in order to have non-singular behaviour. 

IV Discussion 

In this paper we have derived sufficient conditions for an Abelian diagonal or- 
thogonally transitive spacetime with spacelike orbits and with a stiff perfect 
fluid as matter content to be geodesically complete. One of the metric func- 
tions appears to be determinant for the regularity of the spacetime. These 
conditions are easy to check and do not mean much restriction on these 
spacetimes. 

This means that non-singular spacetimes are not as scarce as it was 
thought, considering the reduced list of geodesically complete perfect fluid 
cosmologies in the literature. Further work is needed with more generic sym- 
metries and matter contents in order to clarify the issue, since stiff perfect 
fluids are rather peculiar. They may be interpreted as a massless scalar fleld 
and they are the limiting case for which a barotropic perfect fluid with lin- 
ear equation of state satisfles every energy condition. These spacetimes also 
fulflU the generic condition and are causally stable. 

The latter assert is true since they possess a cosmic time, which is the 
coordinate t. This coordinate has a timelike gradient everywhere. Therefore 
j2], these spacetimes satisfy weaker causality conditions. For instance, the 
chronology condition is true for them and no closed causal curves are possible. 

As it was stated in the introduction, the existence of these non-singular 
spacetimes is possible because they do not possess causally trapped sets. 
They obviously do not contradict then the singularity theorems. They just 
fall out of their scope. 

Another interesting point that is worthwhile mentioning is that the regu- 
larity theorems appear to encourage a growing K for large values of \t\. This 
seems to support a conjecture that states that the spatial average value of 
the pressure in non-singular spacetimes is zero [TT], since p decreases with 
large K according to (f24j) . 



In our regular spacetimes, constant t sheets are Cauchy hypersurfaces 
and we may write the whole system of equations as an initial value problem 




(33) 
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for U, K and H for any constant t. Without breaking the generahty of the 
result, we may focus on t = 0. The initial value problem can be stated as 



Utt -Urr- — ^0, (34a) 
r 

CTtt -CTrr ^ = 0, (34b) 

r 

Kt^Ut + 2r {UtUr + (Tt(Tr) , (34c) 

C/(0,r) = /(r), Ut{0,r) = g{r), (34d) 
a{0,r)^Uir), ^7^(0, r) = ^,(r), (34e) 
KiO,r)^hir), (34f) 

and the remaining equation in the system, 

Kr^Ur + r{U^ + + al + a^) , (35) 
is used to complete the initial data, 

h{r) = C/(0,0)+ / dr'Kr{Q,r') 
Jo 

= /(r)+ r dr'r'{g{ry + f{ry + g,{r'y + f,{rf}. (36) 
Jo 

In order to know the pressure on the hypersurface t — we just have to 
prescribe the initial data, 

p(0,r) = e-2'^«{^.(r)2-/;(r)2}. (37) 

We show that it must necessarily vanish at infinity if the spacetime is 
causally geodesically complete. 

If the term — f^in the pressure does not tend to zero at infinity, the a 
terms would contribute to h as (if gl + tends to a constant) or greater. 
Unless a negative / overcomes this quadratic term, we would have a pressure 
decreasing as a Gaussian exponential and the average on t — would be 
zero. 

But / cannot beat a quadratic term, because the /' term in the integral 
would mean a positive contribution to h, and we would have again a 
negative exponential. That is, if g^ — does not vanish at infinity, it grows 
much slower than the exponential term decreases and the pressure tends to 
zero. 
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The only possibility we have left then is a positive exponential. This 
means a negative h. If we want / to overcome just the /' term in h, we 
require \f{r)\ < Inr for large values of r, a very narrow strip. 

But we also need to keep under control the a terms in h. They remain 
bounded for large values of r if r'^{g'^ + f^) tends to zero. This means that 
the a term in the pressure decreases faster than r~^. Admitting that h{r) 
might behave as — Inr for large r, the exponential in the pressure would be 
a term, that cannot compensate the a term. 

Therefore, pressure tends to zero for large r on constant time hypersur- 
faces, thereby supporting Senovilla's conjecture in [TT] . 
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